Inelastic electron transport in granular arrays 
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Transport properties of granular systems are governed by Coulomb blockade effects caused by the 
discreteness of the electron charge. We show that, in the limit of vanishing mean level spacing on the 
grains, the low-temperature behavior of Id and 2d arrays is insulating at any inter-grain coupling 
(characterized by a dimensionless conductance g.) In 2d and g ^> 1, there is a sharp Berezinskii- 
Kosterlitz-Thouless crossover to the conducting phase at a certain temperature, Tbkt- These 
results are obtained by applying an instanton analysis to map the conventional 'phase' description 
of granular arrays onto the dual 'charge' representation. 

PACS numbers: 73.23.-b, 73.23.Hk, 71.45.Lr, 71.30.+h 



I. INTRODUCTION 

It has been long appreciated that the low-temperature 
physics of generic disordered metals is characterized by 
a subtle interplay of electron-electron interactions and 
coherent disorder scattering. While both effects are of 
crucial importance, their unified treatment still evades a 
complete theoretical description. It is useful, thus, to ap- 
proach them separately. The limiting case of "coherence 
without interactions" has been studied intensely. It is 
well understood that the coherent multiple scattering off 
impurities leads to Anderson localization: in one and two 
dimensions all states are localized 1 -. While in homoge- 
neously disordered systems this phenomenon has always 
to be taken into account, granular systems admit for a 
parameter regime where the physics is entirely controlled 
by interaction effects. It is the purpose of this paper to 
explore the regime of "interactions without coherence" 
accessible in metallic granular arrays. 

Metallic granular arrays are also of great interest 
in their own right 2 i 3 i 4 i 5 i 6 i 7 i 8 i 9 i 10 i n i 12 i 13 i 14 i 15 . In par- 
ticular, strongly coupled arrays (g 3> 1, where g is 
the dimensionless inter-grain conductance) have become 
a subject of increased theoretical attention in recent 
years&SiiSiiiiiSiii^ii. An isolated grain or quantum dot is 
characterized by three energy scales: the Thouless energy 
£/rh, the charging energy E c , and the mean level spac- 
ing S. The tunneling coupling between the grains adds 
another parameter: the dimensionless conductance g of 
a contact between two neighboring grains. Throughout 
this paper we focus on the regime, where the Thouless 
energy £rh is the largest energy scale. This allows one to 
treat each grain as a zero-dimensional object and disre- 
gard the intra-grain in comparison with the inter-grain 
resistance. The interaction effects are controlled by the 
charging energy E c (in the simplest model E c = e 2 /(2C), 
where C is the self-capacitance of a grain.) In our studies 
it is the next largest energy scale in the system. Finally, 
quantum coherence effects are governed by the energy 



scale(s) ex S. If such a scale is much smaller than all rel- 
evant temperatures, one may treat each grain as having 
a continuous spectrum. This assumption allows one to 
disregard phase coherence. In essence: an electron exit- 
ing a grain is never the same electron that has previously 
entered it. 

The parameter regime specified above justifies the 
"interactions without coherence" program. It is clear 
though, that such a simplification cannot work down to 
the very smallest temperatures. At low enough temper- 
ature, coherent propagation through multiple grains will 
become important and our approximation is bound to 
fail. It is, thus, important to realize that the subject of 
our considerations is a transient (though possibly wide) 
temperature range. In this range, coherence may be 
disregarded while interactions (and inter-grain tunnel- 
ing) are crucially important in determining the electrical 
properties of the array. 

At small inter-grain conductance g <C 1, an electron is 
completely localized within a single grain. Therefore, the 
problem is reduced to the description of classical charges 
moving on the lattice (which is a simple limiting case of 
the considerations given below.) In the simplest case of 
on-site interactions only, there is an energy barrier E c 
impeding the transition of electrons between two neigh- 
boring grains. It is thus natural to expect activation 
behavior of the conductivity, with the activation temper- 
ature T* = E c . 

The present paper is devoted to the more intricate case 
of large inter-grain conductances, g ^> 1. In this case the 
charge may spread over many grains to decrease its charg- 
ing energy. The interplay of interactions and tunneling 
dictates that this spreading involves an (exponentially) 
large, but finite number of grains. As a result, the lowest 
energy excitations of the system are large single-charge 
solitons. The activation energy for creating such an ex- 
tended charge carrier is substantially reduced, leading to 
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the low-temperature conductivity of the form: 
v( T ) = 9 exp 



T* 



(1) 



where T* ld) ~ gE c exp[-g/4] and T* 2d) ~ 5 2 £ c exp[-g]. 
The important consequence of Eq. (JTJ is that one- and 
two-dimensional arrays are insulating at arbitrarily large 
inter-grain conductance, g. This is a pure interaction ef- 
fect; Anderson localization physics is not included in the 
model. Switching off the interactions, one obtains Ohmic 
metallic behavior with a temperature-independent con- 
ductivity. 

The solitons interact with each other up to distances 
comparable to their (exponentially large) radius, even 
if the initial model possesses on-site interactions only. 
Once they start to overlap, Eq. Q is not valid anymore. 
In Id this happens at T ~ T{* ld y where the conductivity 

smoothly crosses over to its high-temperature behavior 9 
07x<fl(T) = g — 2ln(gE c /T). In two dimensions, the soli- 
tons interact logarithmically over a large range of dis- 
tances. This leads to a Berezinskii-Kosterlitz-Thouless 
(BKT) unbinding of soliton-anti-soliton pairsiSiii at the 
temperature 



?bkt = TT 2d) /g < T/* 2 



" (2d) 



(2) 



Around this temperature the conductivity undergoes a 
sharp crossover from the exponentially small value given 
by ct(Tbkt), cf. Eq. (|TJ, up to the high-temperature 
asymptotics, ov 2( «(T) = g — \n{gE c /T). In the model 
with only mutual capacitances between neighboring 
grains the Coulomb interaction is logarithmic at arbi- 
trarily large distances. This results in a true BKT phase 
transition with zero conductivity below the transition 
point. The j < 1 version of the latter model was pre- 
viously considered in Refs. |3ll8|. The introduction of 
on-grain Coulomb interactions transforms the transition 
into a crossover. Interestingly, for g ^> 1 the BKT re- 
mains sharp even for the pure on-grain (self-capacitance 
only) Coulomb interactions. 

Technically we approach the problem from two com- 
plimentary perspectives: the phase and the charge rep- 
resentations. The former is straightforwardly derived 
from the microscopic fermionic models. It is commonly 
employed in the study of both homogenous and gran- 
ular interacting systems. While being effective in the 
high-temperature regime, it becomes increasingly diffi- 
cult to handle at lower temperatures. To treat this lat- 
ter regime, we employ the charge model, introduced pre- 
viously within the context of quantum dot physics 20,21 . 
Our main technical achievement is the proof of equiva- 
lence of these two approaches over a parametrically wide 
range of temperatures. For these temperatures, both 
models may be handled in a controlled way. We thus 
conclude that the charge model, although not directly 
deduced from the microscopic Hamiltonian, is indeed the 
proper description of the low-temperature phase of the 
system. The results mentioned above (as well as others 



discussed below) then follow in an almost straightforward 
manner from the charge description. 

The equivalence of the two approaches is based on 
a very important observation. The charge discreteness 
(crucial in the low-temperature insulating phase) mani- 
fests itself in the phase model through the 27r-periodicity 
of the phase field (the internal space of the field is the 
circle S .) The latter results in the existence of topolog- 
ically distinct stationary-point field configurations, clas- 
sified by the integer winding numbers W\ (where the 
vector index 1 numerates the grains on the lattice.) In 
strongly connected arrays, j » 1, the action cost for 
configurations with non-zero winding numbers (so-called 
Korshunov instantons^) is exponentially large. How- 
ever, one has to take into account Gaussian fluctuations 
around the topologically non-trivial stationary points, 
which yield a factor (gEc/T) 1 ^ for each winding num- 
ber mismatch between neighboring grains. This factor 
suggests that the instanton configurations are increas- 
ingly important at low temperatures 23 . Summation of 
the instanton "gas" along with the corresponding Gaus- 
sian fluctuations and the phase-volume factors results 
exactly in the classical (low-frequency) limit of the d— 
dimensional charge model. Specifically (see below), the 
instanton expansion of the phase model coincides term by 
term with the perturbative expansion in back-scattering 
amplitudes of the charge model. Therefore, we are con- 
vinced that the explicit account for instantons in the 
phase-like models is imperative to restore the charge dis- 
creteness and, thus, to describe the insulating phase. 

One may justifiably worry about the role of non- 
Gaussian fluctuations. The latter are known to become 
large at a low enough temperature T ~ E c e~ dg / 2 , vio- 
lating the validity of the instanton gas picture. Crucially, 
however, (in Id and 2d) the corresponding charge model 
predicts an activation gap which is parametrically larger 
(exponential (in g) in Id and algebraic in 2d) than To. As 
a result, there is a wide range of temperatures, where the 
fluctuations are well under control, while the physics is 
completely dominated by the proliferation of instantons. 
The latter results in the appearance of the unit-charge 
extended solitons as low-energy charged excitations and, 
thus, in activation insulating behavior, Eq. Q). In 3d, 
proliferation of instantons and the onset of strong non- 
Gaussian fluctuations, resp., take place at comparable 
temperatures. As a result the phase-charge equivalence 
cannot be reliably established. It seems plausible, how- 
ever, that the instanton gas — and thus the correspond- 
ing charge representation — provide a qualitatively cor- 
rect description of the 3d insulator as well. 

This paper is an extension of two previous shorter 
publications (ltHl} . Its intent is two-fold. Firstly, 
we present some new results. In particular, we extend 
calculations beyond the tunneling limit, accounting for 
arbitrary transmission amplitudes between neighboring 
grains. Furthermore, in addition to an evaluation of 
the transport properties, we discuss the behavior of the 
single-particle density of states (DoS). Secondly, we bring 
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out the philosophy of our approach and expose extensive 
technical details of the calculations. Our main message 
is that charge quantization is crucial in describing the 
low-temperature physics of the array — and, therefore, 
a description in terms of charge degrees of freedom is 
appropriate. As mentioned above, this description is ob- 
tained by accounting for topologically non-trivial field 
configurations in the phase picture. This goes beyond 
the commonly used perturbative treatment of the phase 
model. In Id and 2d arrays, the latter completely misses 
the appearance of a new temperature scale T* marking 
the crossover to insulating behavior. 

The paper is organized as follows: in Sec. [HI we intro- 
duce the phase and charge models. Before coming to the 
main part, namely quantum dot arrays, in Sec. IIIII we 
discuss the physics of a single dot connected to two leads. 
Sec. llVl discusses one-dimensional arrays whereas Sec.lVI 
contains the two-dimensional arrays. The conclusion and 
open questions are discussed in Sec. I VII 



II. PHASE AND CHARGE REPRESENTATIONS 

In this section, we introduce two effective models used 
to describe d-dimensional quantum dot arrays. As men- 
tioned in the introduction, the two descriptions are opti- 
mally adjusted to the nominally metallic and the nearly 
insulating regime, respectively. The application of the 
two models to the computation of observables, and the 
mapping of one onto the other will be discussed in later 
sections. 

Widely used in the literature is the so-called 
Ambegoakar-Eckern-Schon (AES) modelSiiS — a de- 
scription of arrays in terms of phase fields. In the limit of 
vanishing level spacing, this model may be derived start- 
ing from a microscopic description in terms of electronic 
degrees of freedom. The model is presented in Sec. Ill Al 
and its derivation is reviewed in App. IA II 

While the AES approach provides an efficient descrip- 
tion of the high-temperature regime, it is untractable 
in the low-temperature regime where interaction effects 
become significant. Rather, at low temperatures, an al- 
ternative description in terms of charge degrees of free- 
dom is more appropriate. This latter formulation may 
be derived from a phenomenological model introduced by 
FlensbergSS and Matvee^iSi. We review the derivation in 

sec, mm 

The equivalence of the two models — established by 
a mapping between them — will be discussed at later 
stages. 



A. Phase model 

In the regime, where the level spacing of the dot is 
negligible, the dot can be described by a single degree 
of freedom. Starting from a description in terms of elec- 
trons, a phase field <f> is introduced to decouple the inter- 



action on the dot. Subsequently the electronic degrees of 
freedom can be integrated out, yielding an effective the- 
ory in terms of </>. The time-derivative of </> corresponds 
to the voltage V on the dot: V(r) = 4>{t), where r is 
imaginary time. 

Since the AES model is largely standard by now, we 
here restrict ourselves to a brief discussion of its main 
elements. (For an outline of its derivation, see Add. I A ll ) 
The phase action S consists of two terms, S — S c + 
St, describing the charging interaction on the grains and 
the tunneling between neighboring grains, respectively. 
For the d-dimensional array geometry, the charging term 
reads 



(3) 



where 1 is a d-dimensional index, denoting the posi- 
tion of the grain. Here, E c — e 2 /(2C) is the charg- 
ing energy, where e is the electronic charge and C the 
self-capacitance of the dot. The dimensionless quantity 
q = V g C/e is the background charge on the dot as deter- 
mined by an external gate voltage V g . The phase fields 
4>\ obey the boundary condition 4>i{j3) — 0i(O) = 2irW 
(where W € Z.) 

The tunneling term is given by: 



i 

Its 



$> fc tr [(Ae^'Ae-^"')' 



<M') 



(4) 



where 4>w = <f>i — <t>v ■ The matrix A takes the form A nm — 
<5„ m sign(e„) in Matsubara basis (e„ = 2n(n + 1/2)T). 
Furthermore, the coefficients Kk are related to the tun- 
neling matrix elements T a and the density of states v as 
Kk = -4^=^ Y. a \xvT a \ 2k . (Note that we do not need 
to require the transmission in every channel, a, to be 
small — tunneling can be taken into account to arbitrary 
order—.) 

Having presented an effective action for the phase field 
4>, we proceed by discussing its properties. A large inter- 
grain conductance g suppresses dynamical phase fluctu- 
ations: in a conductor, voltage fluctuations are small. 
As a first step, one may, thus, expand the action up to 
second order in <f>. The quadratic tunneling action reads 



gT_ 

47T ^ 

<U'},« 



\<hv, 



(5) 



Here, the dimensionless conductance of the contacts, g, 
is given by24 



9 



^2k 2 K k = 22%,, 



(6) 



where T a = 4tt 2 ^ 2 |T q | 2 /(I + 7rV|T Q | 2 ) 2 is the transmis- 
sion probability in the channel a. The action (jSJ) de- 
scribes Ohmic dissipation. Evaluating transport proper- 
ties of the array within this approximation, one obtains 
the classical Kirchhoff laws. 
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Going beyond the quadratic approximation^, however, 
one finds that the conductance is renormalizcd to smaller 
values upon lowering the temperature. Taking into ac- 
count the terms quartic in (f>, one obtains the renormal- 
ized inter-grain conductance, 25 

0-S-^£T a (l-T a )ln^. ( ? ) 

" a 

In the tunneling limit (T a <C 1), Eq. Q reduces to^ 
g(T) = g-(2/d)\n(gE c /T). 

Eq. (jJJ) states that in all dimensions interactions gen- 
erate logarithmic corrections to the inter-grain conduc- 
tance g. This result holds as long as the corrections are 
small. Perturbation theory breaks down at the temper- 
ature where the corrections become of order of the bare 
conductance or, in other words, the renormalized con- 
ductance g(T) reaches values of order 1. This defines 
a temperature scale T ~ E c e~ dg l 2 . The temperature 
range below this scale is beyond the applicability of the 
perturbative treatment. 

However, there is more to be extracted from the AES 
model even at T » T . Let us return to the full action 
Eqs. © and 10}. The phase field (j>\ is a periodic vari- 
able </>i(/3) — </>i(0) = 2nW. Consequently, the conjugate 
variable — which is charge — is quantized. By using 
the perturbative expansion in (f>\ around the minimum 
<f>\ = 0, all information about periodicity and, therefore, 
about charge quantization is lost. Although for g 1 
phase fluctuations are heavily suppressed, there are <p\- 
configurations that explore this periodicity and, thus, in- 
corporate the manifestations of charge quantization. As 
we show in this paper, these phase instanton&& provide 
the key to access the insulating phase of arrays. 

Before discussing the instanton physics in single quan- 
tum dots as well as one- and two-dimensional arrays, we 
introduce the aforementioned alternative model to de- 
scribe the system. 

B. Charge model 

In this section, we introduce a phenomenological model 
describing the system in terms of charge degrees of free- 
dom, i.e., those degrees of freedom that become ap- 
proximately conserved in near insulating regimes and, 
therefore, are optimally suited to describe the low- 
temperature physics of the system. 

Let us start by considering a point contact between a 
quantum dot and a metallic reservoir. Due to size quan- 
tization effects, no more than a few transverse modes 
a = 1, . . . , N are permitted to transport charge across 
the contact. Each of these modes may be thought of 
as a one-dimensional electron liquid. For simplicity, we 
focus on the case of just one propagating mode N = 1 
throughout. The generalization to multi-mode contacts 
— essential in order to describe the case of large dimen- 
sionless inter-grain conductance g 1 — is discussed in 
appendix IB II 




lead dot 



FIG. 1: (Color online) Schematic representation of the charge 
model. The line in the middle indicates a clean (infinite) Id 
channel with a single impurity, having backscattering ampli- 
tude r. 



Bosonizing the one-dimensional electron liquid in the 
conventional waj^ the system is described in terms of 
a bosonic field 8(t, £). The gradient of this field, d^9(r, £), 
defines the local electron density, i.e., the electron num- 
ber on the dot may be written as M = f^°d£, d^9(r, £) = 
—9(t, 0). This implies that the Coulomb energy takes the 
simple form {eM) 2 /{2C) = E c 6 2 (t,0). Finally, account- 
ing for backscattering by introducing a point scatterer of 
reflection amplitude r at coordinate £ = (cf. Fig. QJ, 
the imaginary-time action of the bosonic field reads 

S[9(t,z)} = j dr[ J d£ [(d T 9) 2 + (d^) 2 ] + 

-oo 

Dr i 
+E c 9 2 (t,0) cos[27t0(t,O)U, (8) 

7T J 

where D is the electronic bandwidth. Integrating over 
the fields #(£ ^ 0) — which is possible because their 
action is quadratic — we obtain 

P 

m 

(9) 

as the effective action of a single remaining degree of 
freedom 0(t) = 6(t,0). Here we have introduced the 
Matsubara representation 9 m — dr 9(r)e^ luJmT , where 
Lu, m = 2-kTvtl. The dissipative term, ^ujm^^, appears 
as a consequence of the assumption that the mean level 
spacing is the smallest energy scale in the model, 5 — * 0. 
It is generated by integrating out the continuum spec- 
trum of the degrees of freedom on the dot. 

The above expression Eq. © can be easily general- 
ized to the array geometry. A field 6*j,i is assigned to 
each contact, where the d-component index 1 denotes 
its position within the array and i = 1 . . . d labels its 
direction. In this notation, the instantaneous electron 
density on the grain 1 is given by the lattice divergence, 
M = 5^j(#i,i+ei — = V • 9\ (where is a unit vector 
in i-direction and the vector notation 9\ is introduced.) 
The generalization of Eq. © reads: 
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where D is again the bandwidth. As in the single dot ex- 
pression Eq. @ , the first term in the action (|10f> describes 
the dissipative dynamics originating from integrating out 
degrees of freedom within the grains, the second term is 
responsible for the interaction effects, i.e., charging, and 
the third one describes backscattering in the contacts. 
Furthermore, we introduced the external gate voltage, q, 
as an additional control parameter. 

As shown in appendix IB II generalization to the N- 
channel case amounts to replacing the single reflection 
coefficient r by the product n^Li r aj where r a is the 
reflection coefficient in channel a. We define the dimen- 
sionless parameter Go as 



Go = -}]la\r 



_ p -6o/2 



(11) 



Although the single channel expression Eq. (|1C)|> was ob- 
tained for a small reflection coefficient r <C 1 , its multi- 
channel generalization remains valid as long as Go ^ 1 
(i.e., individual reflection coefficients r a may be arbi- 
trary, < r a < 1.) For a tunneling contact T a < 1, 
where T a — 1 — \r a \ 2 is the transmission coefficient in 
channel a, this product can be expressed through the di- 
mensionless conductance g — ^ Q T a within exponential 
accuracy as 



n 



exp 



exp 



In this regime, the dimcnsionlcss parameter, Go — 9- 

The action (jf OJI will be our starting point for explor- 
ing interaction effects in single dots as well as in arrays. 
We will use it as an alternative to Eqs. and @ and 
discuss the connections between the two descriptions in 
the following sections. 



III. SINGLE QUANTUM DOT 

The simplest setup on which the impact of interactions 
on transport through an almost open system can be stud- 
ied is a single quantum dot coupled to two leads^Si. 
Interesting in its own right, the discussion of the quan- 
tum dot will facilitate the development of the formalism 
required to describe arrays. We follow a three-step pro- 
gram: in Sec. If 11 Al the AES phase model is investigated, 
in Sec. If If Bl the alternative charge description is used, 
and in Sec. IIII CI the two procedures are compared. 



A. Phase model 

Consider a single quantum dot coupled to a left (L) and 
right (R) lead. In the limit of a vanishing level spacing 
6 — > 0, the system may be described by the phase action 
Eqs. ® and Q}, 



S = 




(12) 



The perturbative results one may derive form this action 
have been discussed in the previous section. Going be- 
yond this level, we here include topologically non-trivial 
excitations and discuss the resulting charge quantization 
effects. 

In addition to the constant solution = 0, the tun- 
neling part of the action is stationary on the so-called 
Korshunov instanton configurations^. These additional 
saddle point solutions are characterized by their winding 
number W = (4>(f3) — (f)(0)) / (2ir) and can be represented 

^22.27.28.29 



\w\ 

J4>w({z},t) _ J~J 
a=l 



2ltlTT 



-i sign W 



1 



= 2irirT 



(13) 



Here the \W\ complex parameters z = (zi . . . z\w\) are 
subject to the condition \ z a \ < 1. The temporal variation 
of (f> corresponds to a voltage pulse V = i(j> on the dot: the 
parameters 1 — \z\ determine the duration of the voltage 
pulse and argz its instance (z — corresponds to a linear 
phase profile or a constant voltage.) 

In the limit T <C E c , the action is dominated by the 
tunneling term, implying that the Korshunov instantons 
are approximate saddle point configurations of the total 
action. Substituting Eq. (| 1 3|) into the action (|I2|I . one 
finds thatSS 



S[<f> w ({z})]ng \W\-2iriWq, 



(14) 



where the dimensionless conductance Go is given by 
(cf. Eq. ifrTjll Go = J2k K 2fe-i- Apart from a small charg- 
ing contribution, S c [<t> w ({z})] = n 2 {T/E c ) £ a Q ,(l - 

|za| 2 |z a '| 2 )/((l - z a z*,)(l - z*z a >)), the action is inde- 
pendent of z a , i.e., the variables z a are instanton zero 
modes. 

Due to the largeness of the parameter 5o > 1) the 
contribution of a single instanton is exponentially small. 
However, as will be shown in the following, fluctuations 
around the instanton trajectory increase with decreas- 
ing temperature 2 - 8 , i.e., a temperature scale exists below 
which the instanton contributions become important. 
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The partition function Z can be represented as a sum 
over different winding number sectors: 



X -> Z W 2-rriWq 



(15) 



where Zw is the contribution from configurations with 
winding number W . Note that a given total winding 
number W can be obtained by superposition of a se- 
quence of s + W instantons and s anti-instantons. Al- 
though these configurations are not true saddle point so- 
lutions, it can be shown that the interaction between in- 
stantons is weakSI and the ideal (instanton) gas approx- 
imation may be used. Referring for a detailed account of 
the computation of the corresponding fluctuation deter- 
minant 28 to App. IA 21 we here sketch the main steps. 

Starting from the action (|12f> . one expands in small 
fluctuations 5<f> around the instanton configuration 
4>w{{z})- We denote the Gaussian fluctuation contri- 
bution to the action by 

£Sin S t = g{8(j>\F w \84>), 

where the linear operator Fw is specified in the appendix. 
The spectrum of Fw is given by 



0. 



\w\, 



1 < \m\ < \W\ 
\W\ < \m\. 



To find Zw, one has to integrate over the massive modes 
with eigenvalues X^m^ as well as over the zero modes z a . 



The massive mode integration leads to a reduction of the 
instanton action, 



Sir, 



Go\W\ 



(Go 



(16) 



This "renormalization" of the coefficient Go is analogous 
to the renormalization of the conductance in the Ohmic 
model discussed in Sec. Ill Al (see Eq. 0.) In the present 
context it signals that instantons become increasingly im- 
portant at low temperatures. 

Finally, the integration over zero modes obtains a pref- 
actor ~ (g In E c /T) ' w \. Combining all contributions and 
accounting for combinatorial factors we thus obtain 



E 



i 



Zw 



Zo j^(s + \W\)\s\ 



wg ^ e -So ln E o 



T 



T 



2s+\W\ 

■ (17) 



We finally sum over winding numbers W to obtain the 
instanton contribution to the free energy 2 ^, F = —T In Z, 

5F(q, T) = -2ng 1 E c e- Qa In cos(27nj). (18) 

Charge quantization renders the free energy a periodic 
function of the gate voltage q. However, as expected, the 
amplitude is exponentially small in Go- 

Using the same formalism, we may compute the con- 
ductance G(q) through the dot. As shown in App. IA 31 
the phase representation of the Kubo conductance is 
given by 



J 



G(q) 



^Z-'Ye 2 ^ lim -3 



\{8(f>\F w \m)\ 



') 



s w [&4>] 



(19) 



Summing over instanton configurations, one obtains 

= 1^1-^006(27^, (20) 

where E c — g 2 E c exp[— Go] may be interpreted as an ef- 
fective charging energy 2 ^!. As with the free energy, the 
conductance contains a weak gate voltage periodic mod- 
ulation. Notice, however, that there is no zero-mode 
factor ~ In E c /T. Rather, the massive mode integration 
leads to the much stronger divergence E c /T. 

The above approach is valid as long as the correction 
is small, i.e., T ^> cos(27rg)|. At smaller temper- 



atures, the instanton expansion becomes uncontrolled. 
The q-dependence of the corrections to the free energy 
and the conductance is linked to charge quantization. In 
the following, we shall study the same problem within the 
charge description. We will show that both models yield 
identical results, and point out a number of similarities 
and differences between the two approaches. 



B. Charge model 

Applied to a lead-dot-lead setup, the charge action 
Eq. I|1(J|) assumes the form 
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S[6] 



^E ( E AumW^+E, 



'L.m 



Or,- 



q5; 



mfl) 



■i=L,R. 



i=L,R 



Ti \ dr cos(27T0i(r)). 



(21) 



Generalizing to multi-channel contacts, we replace — > 
exp[— Gi,o/2]- As long as Gi,o ^> 1, one may treat the 
cosine-term perturbatively. Expanding to lowest non- 
vanishing order in e~ ex - °/ 2 and e~ eR '°/ 2 , one obtains the 
contribution to the partition function 

= ^ e - e °cos(2^) / drdr' SR (>-(<Mr)-M-'))\ 
Z 2tt 2 J \ I 

where Go = {Gl.o + Gr,o)/2. 

Using that (e lX ) — exp[— 5 (A 2 )] and evaluating the 
correlators (OiOj), we thus obtain 

Z PI 

— ~ E c e- g ° cos(2irq) / dr dr 1 r (22) 

Z J \t-t'\ 

= ^e- 5o ln^cos(27r g ). (23) 

This correction corresponds to the one-instanton correc- 
tion to the phase model. In a similar manner, the in- 
clusion of higher order terms yields the correction to the 
free energySi, 

SF(q,T) = -^^ e -So ln ^ cos( 2^) ; (24) 

where C « 0.577 is the Euler constant. 

The phase-dependent contribution 5F matches the re- 
sult obtained from the phase model, Eq. i|18|) . Again, 
however, the approximations leading to this result are 
limited to temperatures higher than a certain cutoff tem- 
perature. Presently, however, we are in a position to ex- 
plore what happens below that temperature: 

At smaller temperatures, the cosine-potential itself 
provides a mass to the field 9. The relevant scale may 
be extracted by comparing the amplitude of the cosine- 
potential with the effective bandwidth — they become 
comparable at T ~ E c (q) = E c e~ Go \ cos(2?rg)|. This 
temperature scale provides a cut-off for the logarithmic 
temperature-dependence of the free energySi. Thus, at 
temperatures T < E c {q) 1 the correction to the free energy 
saturates at SF(q,E c (q)). 

Corrections to the conductance may be found in a sim- 
ilar way. As a result, one obtainsSi^ 

G = f (l-^e-^cos^g)). (25) 

To exponential accuracy, and using E c ~ E c e~ So , this 
is identical to the result l|2(J|l obtained from the phase 
model. 



C. Comparison 

In the previous two sections, we reviewed Coulomb 
blockade effects for a single quantum dot strongly cou- 
pled to the two leads. 

In Sec. lIII ^1 we started from a microscopically derived 
phase action. We applied an instanton analysis to iden- 
tify the 'effective charging energy' E c = g 2 E c exp[— Go]- 
We also found that perturbation theory is applicable for 
temperatures T > E c \ cos(27rg)|. 

In Sec. IIII Bl the starting point was a phenomenolog- 
ical model for the same system. The derivation of that 
model required the condition of weak backscattering in 
all channels, r„ < 1. To exponential accuracy, a pertur- 
bative expansion in the reflection coefficients obtained 
results equivalent to those of Sec. lIII Al Again perturba- 
tion theory breaks down when the effective amplitude of 
the cosine-potential becomes of order of the bandwidth, 
T ~ E c \ cos(27r<7)|. At lower tempertures, the tempera- 
ture dependence of the free energy saturates. The con- 
ductance, on the other hand, is suppressed and behaves 
aM (T/E c (q)) 2 - 

As we shall demonstrate below, the equivalence of the 
phase and the charge description, respectively, will per- 
tain to array-type geometries. However, as with the 
single quantum dot, the charge model will be prefer- 
able over the phase model when it comes to exploring 
low-temperature properties. Remarkably, for arrays the 
equivalence of the two descriptions may be demonstrated 
in explicit terms (and not just exemplified on specific 
observables as was the case for the single dot.) A sec- 
ond difference to the isolated dot regards the role of the 
cosine-potential of the charge model: While for a dot, 
a strong potential (of the order of the 'bandwidth') is 
required to impede charge fluctuations, in an array ar- 
bitrarily weak periodic potentials suffice to 'pin' spatial 
modulationsiSiii. In fact, the efficiency of even weak po- 
tentials in impeding charge fluctuations will be at the 
root of the localization phenomenon in arrays. 



IV. ID ARRAY 

We next advance to the Id array geometry. As in the 
previous section, we start with the AES phase model 
(|IV A|) and then continue to the charge model IjIV B|) . 
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Phase model 



Our starting point is the generalization of Eqs. and 
Q to the Id array geometry, 



M /jo 



M-l 



= ^EE K ' tr (Ae«^Ae-«*0 ,(27) 



i=i k 



(26) 



(Ignoring the charging term, the action of a single plateau 
with height \W - W h \ = 2 equals that of two \W- W h \ = 
1 plateaus. However, the latter configuration has the 
benefit of a much larger entropy which is why we will 
ignore the unlikely formation of configurations with step 
heights \W — Wb\ > 1 throughout. The action of the 
phase model for a given configuration of winding numbers 
reads 



(28) 



where M is the number of grains, and 5<f>j — <j>j+\ — 4>j- 
As in the single dot case, we shall focus on instanton 
phase configurations. Consider a configuration, where 
the phase on grain jo has an instanton, i.e., a winding 
number Wj — 1, whereas Wj = everywhere else. The 
action of this configuration is given by St — 2{Go/2) from 
the tunneling terms jo — > jo i 1 and S c = ir 2 T/E c from 
the charging term. Importantly, the fact that the tun- 
neling action depends only on phase differences implies 
the existence of other instanton configurations of similar 
action: for example, the action of a sequence of W = 1 
instantons on the grains {jo, jo + L} will again cost tun- 
neling action St = 2(Qq/2) (from the tunneling terms 
jo - 1 ->■ jo and j + L -> j +L+ 1.) Only the charging 
action depends in an extensive manner on the length L of 
the plateau: S C (L) = Ln 2 T/E c . However, n 2 T/E c < 1 
whereas Go 3" 1- 



W 



■ m 



FIG. 2: (Color online) A typical configuration of the instan- 
ton winding numbers along the Id array. The edges of the 
plateaus play the role of fictitious charges in the Coulomb gas 
mapping. 

Thus, typical instanton configurations consist of long 
plateaus with a given winding number W that differs 
from the winding number of the background, Wb, by ±1. 



The partition function is obtained by summing over all 
configurations {IT;}, 



*= E 



-S[{W,}} 



x (fluctuation terms). (29) 



One may rearrange the sum by introducing a new set 
of variables oi = Wi — Wi-i € { — 1,0,1}. Boundary 
conditions at the leads require Wo = Wm+i = 0, i.e., 
(Ti obeys the sum rule X^i" 1 °« = 0. The contributions 
to the partition function from configurations {07} can be 
classified according to J2i °i = 2& (where the sum is even 
due to the neutrality condition J^. 07 = 0.) We thus find 



E< 



E 

fe=i 



(30) 



-So/2 



2 k 



E 



sf 1 T,],j> \3-o'Wi"j' 



{{^}IEi^ 2 =2fe} 

where we used that |Wj — — > erf as well as 

I x 



E 



^ 2 = -EE-, 

i 3<i j'>i 



2 ^ 



\3-fWw- 



Since the action for 07 = is zero, it is sufficient to focus 
on non-zero charges <7j = ±1 at positions Xi. Relabel- 
ing the Oj, the summation over different configurations 
is replaced by an integration over the positions xi of the 
charges 07. The corresponding partition function can be 
rewritten as 



E 

k=l 



-S(T)/2 



2 k- 



2k M 



n e ex p 



(fc!) 2 AJ- ^ 

v I i=l Xi = l 



7T 2 T 



2fe 



(-I)^"'|^-^| 



(31) 



i.e., as the partition function of a classical Id neutral 
Coulomb gas with the fugacity e~ Q/yT ^ 2 and effective 
charge irT ' j\[E~ c . Here we have replaced Q — » G{T), an- 
ticipating that the integration over fluctuations around 
the stationary instanton configurations will lead to a 



temperature-dependent renormalization of the conduc- 
tance. As a result of a calculation conceptually similar 
to that for the isolated dot geometry (cf. App. IA 2 bl) we 
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indeed find that the fluctuation determinant is given by 



T = exp 



2M gE c 
'Mr"! ~T~ 



T 

gEc 



2M 
M+l 



= TiH ■ ( 32 ) 



The presence of this factor leads to the renormalization 
Go -> Q{T) =Go- i^\n{gE c /T) implied in Eq. 
Notice that for a long array M — > oo the renormalization 
is more significant than for a single grain, M = 1. 

It is well known (see App. IA 4|l that the Id Coulomb 
gas, Eq. (fTTTf) . is equivalently described by the action of 
the Id discrete sine-Gordon model 



2e- e ( T )/ 2 cos(2^) .(33) 



Introducing the parameter 

7 = 47r 2 . 9 e" e ° /2 < 1, 



(34) 



and using that e" e(T)/2 = {gE c /T)er G °/ 2 this action 
may be reformulated as 



S[d] 



T 



- 6>j) 2 — 7 <-«>S 



7 
2tt 2 



(2tt^)}. (35) 



Before discussing the physics of this model, we will de- 
rive it in an alternative manner, viz. from the charge 
model. This will make the equivalence of the charge and 
the phase description explicit, and elucidate the physical 
meaning of the field 9. 



B. Charge model 

The straightforward generaliztion of the charge action 
Eq. (fTTTf) to a Id array of dots is given by 



'j,m 



qS mfi ) 2 + A^r, 



a 2 \ 



Dr f 

— / dr cos(2n9j) 



(36) 



Thermodynamic properties of the array may be probed 
by differentiation with respect to the gate voltage q. Not- 
ing that q couples only to the static sector of the field 
9jfi = J dr 9j (r) , one may try integrate out all non-zero 
Matsubara components 9j^ n ^Q at an early stage of the 
analysis. To this end let us write 

e,j(T)= 9^ + 59^). (37) 

In the quadratic part of the action 9q and 69 decou- 
ple. Denoting the average over the quadratic (W-action 
by (...), we approximate the functional integral over 
the anharmonic part of the action by (exp (J cos 9)) ~ 
exp J (cos 9), i.e., 

/„„„fo~a\\ ( i2n6n I i2ivSB\ . -%2ixB I -i1-n&6\\ 

(cos(27T0)) m ^ o -* g \ e \ e ) +e \ e )) 
= cos(2t76» ) exp [~2n 2 (59 2 )] . (38) 

Performing the integral over 69, we find 

t — i in 
(» 2 ) = -V V - — - — — j r = -Lln-^— ,(39) 

where E(p) — AE c sva 2 {p/2) is the lattice dispersion. 
Importantly, Eq. (|39|l does not contain temperature- 
dependent infra-red singularities. This provides the a 
posteriori justification of the above integration proce- 
dure. Substituting the <5#-averaged cosine operator back 



into the action of the static field we obtain 

Sci = yE {(Oj+ifi-ho-ct) 2 - 75^ cos(277^,o)},(40) 

where 7 = 2ire c r. As shown in appendix IB II general- 
ization to the multi-channel case amounts to the substi- 
tution 7 = 277e c n Q r a ~ e' Go/2 . 

To exponential accuracy, this result is equivalent to the 
action Eq. (|35|l obtained from the Coulomb gas represen- 
tation of the phase model. 

C. Thermodynamics of Id arrays 

In the previous sections, we have shown that the long- 
range physics of the granular array is effectively described 
by a classical model with free energy F = TS C \, where 

F(q) = E c J2{(d 1+ i - 0, - qf - ^ cos(2^)},(41) 

j 

and the field index '0' has been dropped for convenience. 

This model is known as the discrete sine-Gordon model 
or, in the context of the adsorption of atoms on a peri- 
odic substrate, as the Frenkel-Kontorova models. [The 
Frenkel-Kontorova model describes a harmonic elastic 
chain of "atoms" with stiffness E c , placed on top of a pe- 
riodic "substrate" potential of amplitude jE c /(2ir 2 ). See 
Fig. [3]] In the following, we will summarize the physical 
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properties of this model, and translate to the context of 
the Id quantum dot array. 




b) Id array 



FIG. 3: (Color online) Chain of particles (each represent- 
ing one quantum dot) subject to a shallow cosine-potential. 
While an isolated particle would experience no more than 
a weak back-scattering effect, a particle chain is subject to 
strong pinning. This leads to the appearance of a new (pin- 
ning) temperature scale T* 3> To which does not exist for 
single dots. 

In the absence of a gate voltage (q = 0), the ground 
state of the chain is described by 6j — for all j. This 
field describes a zero charge configuration. Charge exci- 
tations correspond to non-vanishing solutions extremiz- 
ing the free energy F. Varying the action one finds that 
these configurations are given by 

2 

6j = - arctan [exp - jo)}] • (42) 

Eq. (|42|) describes a solitary excitation centered at coor- 
dinate jo and extended over a scale £ s = VV7 ^> 1- The 
total charge carried by this excitation is quantized and 
given by Wtot = X^(#j+i -0j) = M +1-&1 = 1< (Due to 
the large inter-grain tunneling,) it is spread over a large 
number -y -1 / 2 ~ e e °/ 4 of grains. 

Substituting this solution into the free energy l|41|) , one 
obtains the soliton energy 

T* = E c ^j » T . (43) 

At finite temperatures, the system will host a gas of 
thermally excited solitons and anti-solitons with density 
n s ~ exp[— T*/T]. Due to the absence of gapless charge 
carriers in the system, transport will exhibit activation 
behavior, as we will show below. 

At finite gate voltage, the uniform configuration 8j = 
will acquire the finite energy 5F e \(q) = ME c q 2 (the elas- 
tic term in the action.) By contrast, the configuration 
8j = jq which minimizes the elastic term has the en- 
ergy 5Fp 0t (q) = —ME c "f/(2ir 2 ) (the cosine-potential.) 
The smaller of the two determines the ground state of 
the system. In the language of the Frenkel-Kontorova 
model, the "incommensurability parameter" q represents 
the periodicity mismatch between the chain and the sub- 
strate. For small values of q the system will find it fa- 
vorable to retain a commensurate state, i.e., the chain 
will stretch a little so as to still benefit from an op- 
timal coupling to the substrate. Thus, the system re- 
mains in a commensurate phase 9j — for gate voltages 
q < q* = \J ~j j (2-k 1 ) . At q* , where the two energies 



become equal (5F e i(q*) = 5F pot (q*)), a commensurate- 
incommensurate transition takes place — at this gate 
voltage solitons are created at no cost. For the average 
number of electrons per grain, Af(q) = q — d q F/ (2ME C ), 
one thus expects: N (q) = for \q\ < q* (insulator) and 
M(q) — > q for \q\ > q* (metal). 




FIG. 4: (Color online) Free energy of the Id array as function 
of the common background charge, q. Only two insulating 
plateaus at q — and q = 1/2 are shown. 

A more thorough discussion of the system (cf. Ref. j33|) 
shows that insulating 'plateaus' along with superimposed 
solitary excitations form not only around q = 0, but 
also around other rational values of q. However, both 
the width of these plateaus and the corresponding ac- 
tivation energies decrease for higher rational fractions. 
Among the low-lying rationals, q = 1/2 plays a particu- 
larly interesting role. Indeed, for a single grain, q = 1/2 
represents the charge degeneracy point, where the sys- 
tem is in a conducting state (Coulomb blockade peak). 
Unexpectedly, the array exhibits a very different behav- 
ior. Using the language of the Frenkel-Kontorova model, 
q = ±1/2 is special in that the atoms of the unperturbed 
chain alternatingly find themselves in minima/maxima of 
the substrate potential. Under these conditions, energy 
can be gained by a 'Peierls-distortion', doubling the pe- 
riod, i.e., a density modulation with amplitude 59 j ~ 7. 
This configuration does not respond to small variations 
in q — corresponding to insulating behavior. However, 
the width of the insulating plateau Aqi/2 ^ 7 is much 
smaller than the width of the plateau Aqo = 2q* ~ ^77 at 
q = 0. The dependence of the free energy on q (including 
the two plateaus) is shown in Fig. 0] 

D. DC transport 

In order to discuss the DC conductivity of the array, 
one needs to restore the low-frequency, uj <C T, dynam- 
ics of the classical charge model, Eq. (|41ll . In principle, 
this may be done by keeping the dissipative term in the 
action. It turns out that in the multi-channel case (cf. 
Add. IB l)l . and for strong backscattering, the coefficient 
of the dissipative term reads as 'Kg~ 1 \ui m \9 2 m . While this 
seems to satisfactorily describe the dynamics of Ohmic 
dissipation, the main drawback of the imaginary time 
approach is that it involves cumbersome analytical con- 
tinuations u) m — > uj + iO. To avoid this complication, it is 
convenient to pass to the Keldysh representation. In fact, 
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it turns out to be sufficient to focus on the semi-classical 
limit of the Keldysh formalism, i.e., a limit physically 
equivalent to a description of the system in terms of a 
classical Langevin equation^. Referring for a more de- 
tailed discussion of these connections to Add. IB 2l we here 
take the adequacy of this formulation for granted and de- 
scribe the system in terms of the Langevin equation right 
away: 

To start with, consider the equations of motion of the 
static model, Eq. gTJ): d (§S0j) - ^ sin(27r0,) = 0. 
Since fjdOj = Vj is the voltage on grain j, this equa- 
tion simply expresses the fact that in the absence of 
charge quantization, 7^0, all grains are equipotential: 
dVj = Vj+i — Vj = 0. We now ask how these equations 
have to be modified if currents are allowed to flow. The 
minimal classical description of dissipative current flow 
between grains j + 1 and j is provided by the current- 
voltage relation Vj+\ — V 3 ■ = Rlj, where R — 2nh/(e 2 g) 



is the contact resistance, and L 



ed t 0j the current flow- 



ing between grains j and j + 1. Restoring the 7-term, 
we are led to the phenomenological generalization of the 
equation above, 



that 9\ has been introduced to describe a change in the 
shape of the soliton in response to the field. This needs to 
be distinguished from the temporal change in its position 
(which has been accounted for by the introduction of the 
as yet undetermined shift vt.) To discriminate between 
these two effects, we require that the linear equation de- 
termining 9\ be an equation in a function space orthogo- 
nal to the zero mode function, 89, describing translations 
of the soliton, 9{x + 5x) — 9(x) ~ Sxd9. In particular, v 
should be determined in such a way that the r.h.s. of the 
equation be orthogonal to 89, 



dx ( -39 
9 



E 



sin 2 (7r6 



= 0. 



(47) 



Using Eq. Q42J1 for 9, this requirement leads to the ex- 
pected result v ~ gE. This in turn implies that the DC 
conductivity is given by a = gn s (T). Using that for low 
temperatures, T < T* , the soliton density shows acti- 
vation behavior, n s (T) ~ exp[— T*/T], we arrive at the 
result Eq. 0. 



- dtO - E c 
9 



d 2 9 - — sin(27r6 
2tt v 



-~E + t(t), (44) 



where we have introduced a Gaussian noise term, £(i), 
with the correlator 



2vrT 



5(t-t')Sj 



(45) 



to satisfy the fluctuation-dissipation theorem, and an ex- 
ternal electric field, fi, as a driving source of current. 
Formally, Eq. l)44f) represents a Langevin equation for 
the classical degree of freedom, 9. 

Our goal is to calculate the current, /, driven by a 
weak uniform field, E. As we saw above, charge trans- 
port in the present model is by solitary excitations. As 
an ansatz for the current we thus use / = en s v, where n s 
is the soliton concentration and v is their effective drift 
velocity. While the soliton concentration, n s , has been 
discussed above, the drift velocity still needs to be deter- 
mined. To this end, we temporarily ignore the noise term 
and seek for propagating solutions of Eq. I|44|). Assum- 
ing the external field to be weak, we consider the ansatz 
9(j,t) = 9(j — vt)+9i(j~vt) + (, where 9 is the static soli- 
ton, 9\ ~ E a small distortion of the soliton shape due to 
the presence of the external field, and the constant £ ac- 
counts for the weak shift of the minimum of the periodic 
potential by the field: -E c 7sin(27r£) = ttE. Linearizing 
Eq. I|44|). we find that 9\ satisfies the equation 



E. Density of states 

As another quantity of interest we discuss the (tun- 
neling) density of states. At large energies, the density 
of states can be conveniently described within the phase 
model. Within the framework of that model, the DoS 



tr [G(iu m )] 



is represented as 



sin ttTt 



*£W>,M-<Mo)) 



.(48) 



At large temperatures, the 'Debye-Waller factor' 
( e <E,(*,(T)-*«(0)^ ^exp[-i£ g (K(T)-<M0)| 2 )] may 
be computed from the quadratic approximation to the 
phase action. This leads to 



1 



(\Mr)-M0)\ 2 ) = 4T 



COS U). m T 



q , m 



MU 2 



E a ■ 



= 2, 



E r T 



9 



(49) 



Integrating over r and performing the analytic continu- 
ation we thus obtain the result 



1/(6) 



vq exp 



) 



(50) 



Ect^Ox = -88-— sin 2 (7r 



9 



(46) 



where = d 2 — 7Cos(27r(9). Importantly, it is not our 
prime objective to identify the actual shape of the soliton 
(#1); Rather, we whish to compute its speed of propaga- 
tion, v. An equation for v may be obtained by noting 



Notice that the conductivity and the DoS, respectively, 
are governed by different energy scales. The DoS be- 
comes exponentially small at energies e ~ E c /g ^> T*. 
To understand what happens below that energy, notice 
that the minimal energy required to add a charge to the 
array is the excitation energy of a soliton. Therefore, at 
zero temperature, the DoS vanishes at energies smaller 
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than T*; i.e. v{e < T*) = 0. The charge quantization 
leads to a hard gap in the DoS. 

This concludes our discussion of the one-dimensional 
array. We have found that the proliferation of instantons 
at low temperatures drives the system into an insulat- 
ing phase where transport is activated and conducted by 
the charge solitons. The temperature at which activation 
behavior sets in is exponentially small in the dimension- 
less conductance. Yet it is parametrically larger than the 
scale To where the perturbative corrections become large: 
T*/T ~ e 5o/4 > 1. The tunneling DoS is significantly 
suppressed at even higher energies and displays a hard 
gap at the soliton energy. 



As a result of a straightforward generalization of the 
one-dimensional fluctuation determinant discussed in 
App. IA2bl we indeed find (cf. App. IA2c|l that the de- 
terminant due to fluctuations around the stationary con- 
figuration is given by 

* - -H^Hm)' ■ ,52) 

We account for the presence of this factor by renormal- 
ization Go -> G(T) =Q - (2/d) \n(gE c /T), where d = 2. 
(In fact, it is straightforward to verify that this result 
holds for systems of arbitrary dimensionality.) 



2D ARRAYS 



B. Charge model 



We next extend our discussion to two-dimensional 
quantum dot arrays. Our strategy will parallel that of 
the previous sections: Starting from the phase represen- 
tation we will establish a connection to the complemen- 
tary charge representation, and then discuss the physics 
of the system in terms of the latter. 



A. Phase model 

Again we start from the action and (0J, where the 
lattice summation now extends over a two-dimensional 
regular array. Anticipating the importance of 'winding 
numbers', we begin with an identification of instanton 
solutions right away. As depicted in Fig. \5\ instanton 
configurations in the 2d geometry assume the form of "is- 
lands", i.e., regions of a certain winding number Wisiand 
surrounded by a background with a different winding 
number W\>- For entropic reasons, configurations with 
Wb = and Wi s i an d = ±1 will dominate. The action of 
one such island of area A and circumference L is given by 



A+^-L. For an island differing by a generic 
winding number from its background, this generalizes to 



7T 2 T 



sland 



-AW 



g{T) 



L\W\ 



(51) 



where we anticipated that the inclusion of fluctuations 
will again manifest itself in a renormalization of the con- 
ductance, Go —> G{T). 



8 v, i+e. 




6v, i 



FIG. 5: (Color online) The "islands" with fixed values of the 
winding number in a 2d array. The area A — 1,2,... is the 
number of grains inside an island, while the circumference 
L = 4, 6, 8, . . . is the number of links crossed by its boundary. 



Starting from Eq. IjlOjl (where the sum now extends 
over a 2d array), our strategy will be to successively in- 
tegrate over high-frequency fluctuations to derive an ef- 
fective low-energy action. (Unlike in the Id case where 
all dynamical fluctuations could be integrated out in one 
sweep.) Consider, thus, fluctuations of 9(t) in a win- 
dow of (Matsubara) energies between the bandwidth D 
and D < D. Integration over these fluctuations will lead 
to a renormalization of the cosine-potential. As long as 
the renormalized backscattering amplitude is small, func- 
tional averages over high frequency fluctuations can be 
taken using the quadratic action. 

Due to the vector nature of 9, there are two types of 
modes contributing to the fluctuations. Using the repre- 
sentation 9 — V% + V x 77, where x an d V are scalar fields 
and V and V x lattice variants of gradient and curl, re- 
spectively, the charging action takes the fy-independent 
form S c [x, v] = E c Xa(^Xi) 2 /T ■ The Gaussian average 



1 

2^ 



E r E 

P D<\u m \<D 



E(p) +n\u; m \ 



thus splits into two contributions, where only one (the 
^-contribution) couples to the lattice dispersion E(p) = 



AE c (sm 



2 p x 



,2 Py 



For frequencies D < E c , we find 



4tt 2 



In- 



D 2 

ErD 1 



(53) 



i.e., the x -m ode is effectively frozen out due to its cou- 
pling to the charging action. However, (and unlike in the 
ld-case) there is one "massless" mode whose fluctuation 
amplitude depends on the effective bandwidth D. Using 
this result to renormalize the coefficient of the cosine- 
potential we obtain 



£) e -So/2 



E r De~ Go/2 . 



(54) 



This integration procedure can be iterated until the 
renormalized backscattering amplitude is of order of the 
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bandwidth, i.e., D ~ E c e~ s< > = To, corresponding to 
an effective conductance G c g — For temperatures 

T > T , (i.e., for a Matsubara frequency spacing ~ T 
larger than the cutoff energy) one may approximate the 
action by the zero Matsubara (classical) contribution 



S c i[&\ 



El 

T 



l{T) 
2ir 2 



J2 cos(27r0 M )}, (55) 



i=x,y 



where j(T) = 2-k^/T] E c e- g °/ 2 . 

At smaller temperatures T < T , all modes become 
massive due to the cosine-term. This leads to a satura- 
tion of the coefficient 7(T), 



\Vt7e~c 



e -g /2 

; e Sa/2 



T>T , 
e -9 T <T Q . 



(56) 



We shall show now that the perturbative expansion in 
powers of 7(T) of the classical model, Eq. (|53)l . leads 
to the same island scenario discussed above within the 
phase model. 

For bookkeeping purposes, we label the coefficients 
7a(T) by the index of the bond a = (i = x,y) 

to which they belong. The partition function can then 
be written as 

oo 

Z = Y / Z n (E c /Tr II 7-(T), 

n— —3 

0=1.. .») 

where Z n is a product of n cosine-terms averaged over the 
quadratic action S c [9] = E c ' 9\) 2 /T. Two types of 

contributions to this expansion can be distinguished: a) 
terms containing higher powers of the cosine taken at the 
same link and b) terms involving different links. One may 
show that the first class of terms describes perturbative 
corrections to the conductance of a single contact (equiv- 
alent to those obtained in Rcf. within the framework 
of the phase model.) We here focus on the phenomenon 
of "island formation" , as described by the second family, 
b). 

The expansion of Z n generates expressions of the form 

exp[ 2 7r*(El=i(±)^x,L fal +E" 1) =,+x( ± )0i/.O]> where 3i 
labels contacts in i-direction. To perform the ^-average, 
we again use the representation fl = Vx + Vx?). While 
the fields \ are defined on the lattice, the fields r\ are 
defined on the reciprocal lattice (see Fig. 0) Impor- 
tantly, the (now classical) rotational field ij is strictly 
massless. This means that contributions to the expan- 
sion containing an uncompensated 77-amplitude in the 
exponent average to zero. The surviving terms obey 
the condition J2 jx HVh„+e y ~ VhJ + Ej y T(r)i h +e x - 
i]i jy ) — 0. This corresponds to the island structure: 
the lowest-order non-local term is proportional to 7 4 = 
7z,i7x,l+e a! 7j/,i73/,l+e ! , — involving all the four links sur- 
rounding grain 1. 

Every island is weighted by a factor (E c j(T)/T) L , 
where L is the order in perturbation theory required 
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FIG. 6: (Color online) The rotational field rj. The lattice curl 
is eiven as (V x r?)i = ( v '+"y Vl ) . Goine around one srain, 

v " V -»71+e x +»7l/ 

the //-fields coming from the 4 links cancel — thus, yielding 
a non- vanishing contribution to the expansion in 7. 



for its formation. For an island with winding number 
I IT I = 1, L is the lattice circumference. Islands with 
\W\ > 1 are surrounded by a chain of 7-amplitudes |IT| 
times; in this case, L = \ W\x circumference. Finally, the 
averaging over the massive x -ne lds results in a factor 
exp[—n 2 TAW 2 /E c ], where A is the island area. 

Summarizing, every island carries a factor Pyy(A, L) = 

{e-^ T ' E *y W {E cl {T)/{2u 2 T)) Lm . Using Eq. ®, 

one finds that in the high-temperature regime, T > 
T , E c7 (T)/(2tt 2 T) ~ e- e ( T )/ 2 . Thus, P W (A,L) = 
exp [-tt 2 {T/E c )AW 2 - \Q{T)L\W\], which is in per- 
fect agreement with the prediction of the phase model, 
cf. Eq. H51[). At lower temperature, non-linear fluctu- 
ation corrections in the phase model diverge^ and the 
non-interacting instanton treatment runs out of validity. 
However, having established the equivalence of the phase 
and the charge model at T > To, one may proceed with 
the analysis of the latter even at smaller temperatures. 



C. Solitons and the BKT crossover 

In order to extract the low-temperature behavior of 
the array, we investigate the properties of the classical 
model (|55|l . The lowest energy configuration of the ac- 
tion is given by the homogeneous solution 8 = (mod 1) 
everywhere. In order to minimize the cosine-potential, 
localized excitations must have integer far away from 
the core. Since the total charge of such a localized excita- 
tion is ej (d 2 l) V • 9 = ej ds- 9, where the line integral on 
the r.h.s. is calculated over a distant contour enclosing 
the excitation, the charge of the excitation is quantized 
in integer multiples of e. Excitations of lowest energy 
have charge ±e. They consist of a large (i.e., spread 



out over 



1/7 



^>> 1 grains — see below) local- 



ized 2d soliton, connected to a Id string of links with 
9i = 1. The other end of the string may either go to 
the system boundary, or terminate in an anti-soliton of 
opposite charge. The soliton solution centered at 1 = 
can be written in the form 9\ = 1 — i?(l) for the links 
along the string and 9\ = §(X) everywhere else, where 
|i?(|l| — > 00) I — * 0. Minimization of the action, Eq. 15511 . 
with respect to •& yields the saddle point equation for the 
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soliton solution, 



V(V-tf) 



2tt 



sin(2vri? 4 )e 4 = 0. (57) 



i=x,y 



Except for a domain consisting of 0(1) links closest to 
the core of the soliton, $ is small, justifying an expansion 
of the sine-term in the saddle point equation, i.e., V(V • 
d) — 7$ = 0. Its unit-charge solution is 



(58) 



where e; = l/l and K\ is a modified Bessel function. The 
large size of the soliton, £ s = 1/ y/j{T) ^> 1, justifies the 
continuum approximation. 

To obtain the soliton energy T* , we substitute this 
solution back into the action, Eq. 155|) . One finds that 
the energy originates primarily from the cosine-potential 
part and is given by T* = (E c j(T)/(2tt)) ln£ s . The large 
factor ln£ s = — |hi7(T) » 1 is due to the logarithmic 
spreading of the charge density over the wide range of 
distances 1 < I < £ s . Due to this factor T* > T ~ 
E c "f(T ) (while their ratio is exponentially large in g in 
Id, here it is only algebraic.) At larger distances, I > £ s , 
the charge density decays exponentially. 

The charge spreading leads to a long-ranged soliton- 
soliton interaction. The solitons interact logarithmically 
up to a distance £ s beyond which the interaction is ex- 
ponentially screened. Since the density of thermally- 
excited solitons is n s ~ exp[— T*/T], the mean dis- 

— 1/2 

tance between them is l s = n s w exp[T*/(2T)] 
which becomes comparable to ( s at T « T*/(21n£ s ) = 
E c j(T)/ (4ir). This condition is satisfied at temperatures 
about the "freezing" temperature, T ~ To. Thus, at 
T < To, the thermally-excited charges are essentially 
non-interacting, while, at T > To, there is a neutral (in 
average) gas of logarithmically interacting solitons and 
anti-solitons. The soliton core energy yields the fugacity 
/ of the logarithmic gas: In / ~ E C ^(T)/T. 

In the latter regime the partition function of the 
charged degrees of freedom may be written as 



n s ~ exp[— 26-\/Tbkt/(T — Tbkt)], where & is a con- 
stant of order unity, driving the array into the conducting 
phase. 

Notice that the Coulomb interactions in our model 
are strictly on-site (only the self-capacitance, C, is in- 
cluded.) The long range of the soliton-soliton inter- 
actions is due to the fact that in a strongly coupled 
array, Go 3> 1, the charge is spread over a large dis- 
tance £ s ~ exp[£/o/2]. To describe the truly long-range 
Coulomb interactions, one may modify the model by 
including mutual capacitances C between neighboring 
grains. It is straightforward to show that such modi- 
fication alters the range of logarithmic interactions as 
£ s — > £ s a/1 + C'/C, while the charging energy now reads 
E c = e 2 /{2(C + C')). In the limit C -> 0, while C 
remains finite, the interaction range diverges, £ s — > oo, 
i.e., one deals with the true (logarithmic) 2d Coulomb 
interaction (without the self-capacitance no electric field 
lines can leave the system.) In this case, one observes a 
genuine BKT phase transition: £ s — » oo and the density 
of free charges below Tbkt becomes strictly zero. 




FIG. 7: (Color online) Phase diagram of the 2d granular array. 
The thick line indicates the BKT crossover temperature at 
g > l. At smaller g the crossover disappears (in the self- 
capacitance model.) The shaded area indicates the region 
with exponentially small conductivity. 



oo 



71=0 



(d 2 l n ) e 



■ 2ttT 



£ ball* 



-VI 



(59) 

The plus/minus signs in the exponent correspond to 
soliton-soliton and soliton-anti-soliton interactions, re- 
spectively. 

It is well known that the Coulomb gas in 2d de- 
scribed by Eq. I|59|l undergoes the BKT transitioniSiii 
at a critical temperature Tbkt ~ £ ; c7(Tbkt)/(47t). 
For T < Tbkt, the charges are bound in charge- 
anti-charge pairs. In this regime, the finite interaction 
range £ s leads to an exponentially small residual den- 
sity of free charges n s w exp[— T*/T], where T* = 
^bkt hi(£ s 2 ). Above the transition/crossover tempera- 
ture the density of free charges rapidly increases asi& 



The conductivity of the array may be evaluated in the 
same way as in the Id case. With the soliton velocity v ~ 
gE, one finds a ~ gexp[—T*/T] at temperatures below 
Tbkt- Above Tbkt, the conductivity behaves 35 as a ~ 
gexp[— 26-\/Tbkt/(T — Tbkt)], whereas, at even higher 
temperatures, this behavior crosses over to the result of 
the perturbative calculation, a = g — hi(gE c /T). 

The corresponding phase diagram is shown in Fig. 
Unlike previous works^iSi 2 ^ that predicted a zero- 
temperature metal for g > g c ~ 1, we find that the low- 
temperature phase is an insulator for arbitrarily large 
g. The critical temperature, Tbkt(s), however, drops 
sharply at g ~ 1 and, at large g ^> 1, behaves as 
Tbkt ~ E c gexp[-Go}. 
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D. Finite gate voltage 

So far we have restricted ourselves to the case of zero 
gate voltage only. A finite gate voltage induces a contin- 
uous background charge q oc V^ a tc on the grains. In this 
case the charging term in the action Eq. H55[) has to be 
replaced with S^[6; q] = (E c /T) ^(V • 61 - qf . Alter- 
natively one may shift the 9 field by ql to move the q— 
dependence into the pinning term cos(2ir(Qi + qU)). 
Since the grain coordinates U take integer values, the 
model is periodic in g-space with unit periodicity. Leav- 
ing the analysis of random background charges for future 
studies we will restrict ourselves to uniform gate voltages 
q(l) = q = const, throughout. 

For small q, the system is in a particle-hole symmetric 
"neutral" state: as for q = 0, the ground state is given 
by 0\ = 0. At some finite value q — q*, a transition 
towards a charged (with a non-integer average number 
of electrons per dot) and spatially non-uniform ground 
state takes place. To find q* , let us compute the soliton 
energy in the presence of q. Since the g-dependence of 
the Hamiltonian is a pure boundary effect, one immedi- 
ately finds the soliton energy T*(q) = T*(0) - 2qE c . At 
q* = T*(0)/(2E C ) the soliton energy T*(q) vanishes. This 
marks the transition into the charged state: for q > q* , 
solitons are created at no cost. 
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FIG. 8: (Color online) The phase diagram of the 2d array with 
a uniform background charge, q. In the "neutral" (insulating) 
phase, there are (almost) no solitons. Outside the neutral 
regime, charged solitons are generated. Below the melting 
temperature T m , they form an insulating Wigner crystal. 



It is instructive to relate the physics of the array to 
the phenomenon of vortex formation in type II supercon- 
ducting films 36 . Within this analogy^, the gate voltage 
translates to an external magnetic field, H, and the gate 
voltage q* corresponds to the critical magnetic field H c \ 
where vortex formation becomes energetically favorable. 

Above H c i a type II superconductor contains a fi- 
nite density of vortices which at low enough tempera- 
tures form an Abrikosov lattice^ 9 . In a clean film, the 
vortex lattice is free to move, but is easily pinned by 
the system boundaries, the underlying lattice structure 



(as in Josephson junction arrays) or any sort of disor- 
deip2L££. Upon increasing the temperature it will even- 
tually melliiSii ., and above the melting temperature T m 
most of the vortices are unbound. The melting tempera- 
ture at finite H > H c i is smaller than the zero magnetic 
field Berezinskii-Kosterlitz-Thouless but parametrically 
of the same order—. Thus, at T < T m the system is su- 
perconducting while at T > T m vortex motion leads to 
dissipation. 

Translating back to our problem this means that at 
q > q* the solitons form a Wigner crystal once their den- 
sity is sufficiently large such that the interaction is loga- 
rithmic. Only in the narrow interval q* < q < q* + £.7 2 
the system is in the conducting charge liquid state. Upon 
increasing the gate voltage, the Wigner crystal forms and, 
due to lattice pinning, the system is an insulator at tem- 
peratures smaller than the melting temperature. The lat- 
ter is of the order of Tbkt ■ The phase diagram is shown 
in Fig. El Note that while for q < q* charge is carried by 
individual (thermally-activated) solitons, for q > q* the 
mobile charges are lattice defects, whose core energy is 
proportional to the logarithm of the lattice constant of 
the Wigner crystal. 



E. DoS 

As in the Id case, interactions lead to a suppression 
of the tunneling DoS in the vicinity of the Fermi energy. 
Using Eq. I|48|l . one finds that, in 2d, the suppression of 
the DoS is determined by the phase correlator 



1 — COS Ul m T 



2irg 



\n 2 (gE cT ) 



(60) 



The final expression for the DoS, thus, read; 



,9,41 



u(e) ~ vq exp 
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Ang 



In" 



gE c 



(61) 



As in one dimension, the suppression of the DoS becomes 
significant at a scale much larger than T*, namely, at 
e ~ gEffi^" 2 ^ 1 ^ . While this is a purely perturbative 
result, decreasing the energy further non-perturbative 
effects become important. Specifically, the finiteness 
of the soliton energy leads to a hard gap in the zero- 
temperature DoS, 



T=0 



is{e<T*) = 



VI. CONCLUSIONS 



(62) 



We have studied transport properties of inelastic gran- 
ular arrays. We find that charge quantization and the 
localization of a unit charge within a finite area play the 
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crucial role in the low-temperature behavior of these sys- 
tems. This is the case even in systems where due to 
a large bare conductance g ^> 1 charge behaves like a 
fluid at high temperatures. At large g, the elementary 
charged excitations are solitons spreading over (exponen- 
tially) many grains. In contrast, for weakly coupled ar- 
rays (j< 1) charge is quantized on a single grain. The 
excitation energy of the solitons, T* , determines the ac- 
tivation gap for the low-temperature conductivity. In 2d 
arrays, logarithmic interactions between solitons lead to a 
sharp BKT crossover to the high-temperature regime at 
a temperature Tbkt which is parametrically smaller than 
the activation gap: Tbkt = T? 2d Jg -C TT 2d y The most 
straightforward way to access this physics is to study the 
elementary excitations of the classical, dissipative charge 
model. The latter may be viewed as the low-frequency 
limit of the phenomenological charge model, Eq. (|10fl . 
introduced in the context of a single do1iS2i2i. The re- 
duction to the low-frequency (classical) sector proceeds 
through integrating out all the modes with frequencies 
higher than a certain bandwidth (eventually to be under- 
stood as the temperature.) This process renormalizes the 
initial amplitude of the cosine-potential E c e~^°/ 2 , bring- 
ing it down to E c e~ Ga/2 {T/E c ) 1 - 1 / d . (The renormaliza- 
tion is due to the massless rotational modes, which rotate 
charge without compressing it. Such modes are absent 
in Id, thus there is no renormalization for d = 1.) The 
effective bandwidth reaches the amplitude of the cosine- 
potential at a certain energy scale (freezing temperature) 
To. The latter is determined by the condition: 

Go - ~ HEc/To) = . (63) 
a 

It is reasonable to assume (based on the exact solution 
of the single-dot model--) that the renormalization of 
the potential stops at this scale. The physical reason 
behind this saturation is freezing-out of the rotational 
modes due to mass generation (by the cosine itself.) It 
is very important, however, to realize that in Id and 2d 
the soliton energy is larger than the freezing tempera- 
ture: T* ^> T . As a result, the insulating behavior is 
established independently of the validity of the freezing 
assumption (the latter simply allows one to describe the 
insulator at T < T .) 

In this paper we have derived the above mentioned 
phenomenology starting from the complimentary phase- 
model. The latter is directly deduciblc from the micro- 
scopic fermionic Hamiltonian. The mapping between the 
two approaches is achieved by summation over instanton 
configurations (plus Gaussian fluctuations) of the phase- 
model. Curiously, the non-Gaussian fluctuations of the 
phase field around the instantons become strong at the 
same temperature scale To given by Eq. H63|) . As was 
already mentioned above, the insulating behavior in Id 
and 2d sets in at the scale T* , which is parametrically 
larger than To given by Eq. (|63|) . Thus, the instanton 
gas summation is fully justified. In 3d, T^ 3d ^ w To, and 
therefore instantons and non-Gaussian fluctuations be- 



come important at the very same temperature, compli- 
cating the theoretical treatment. Our approach allows 
one, thus, to follow the behavior of the (Id and 2d) sys- 
tem from the metallic phase at high temperatures down 
to the insulating phase at low temperatures. 

However, this luxury comes at the price of several key 
simplifications of the model, (i) The array is assumed 
to be perfectly uniform, i.e., we choose grain capaci- 
tances, C, inter-grain couplings, g, and (most restric- 
tively) background charges, q, to be identical for all the 
grains, (ii) We have assumed a continuous spectrum in 
each grain. Doing this we have disregarded all manifes- 
tations of quantum coherence, such as Anderson local- 
ization, elastic hopping through several grains, etc. As a 
result, our treatment misses the physics at energy scales 
associated with the single-particle level spacing, 5. Thus, 
it provides only transient temperature dependencies - 
and not the ultimate low-temperature conductivity. 

In principle, (random) fluctuations of the parameters 
(C, g, q) may be incorporated in our treatment. They 
will directly translate into static fluctuations of the cor- 
responding parameters of the charge model. If such fluc- 
tuations are smooth and long-range correlated, solitons 
will continue to be the elementary charged excitations of 
the model. The difference is that now the solitons move 
in the presence of a random pinning potential. The prob- 
lem is reduced to the dissipative dynamics of a classical 
gas of interacting solitons and anti-solitons subject to 
pinning. In 2d the very same model is used to describe 
the motion of vo rtices in type II superconductors, see, 
e.g., Refs. |38l42| . In Id the conductivity of the array 
is determined by rare fluctuations causing anomalously 
strong soliton pinning 43 . Strong short-range disorder, 
however, may invalidate the soliton picture. It is not 
known whether arrays with short- and long-range disor- 
der display the same low-temperature behavior. 

Incorporating quantum coherence is yet a more diffi- 
cult task. To proceed in this direction, one has to take 
into account the finite level spacing 5 in each grain. The 
most natural model assumes a random chaotic spectrum 
in each grain with Wigner-Dyson spectral statistics, mu- 
tually uncorrelated between different grains. Such a 
model necessarily includes randomness, even if all the 
other parameters are assumed to be strictly determin- 
istic. Technically, the chaotic spectrum of a grain is 
described by the non-linear cr-model matrix field Q\. 
Since our ultimate goal is to describe the interacting sys- 
tem, either replica^ 4 , or Keldysh 45 variants of the latter 
should be considered. In this case the Qi-field is a ma- 
trix in replica (Keldysh) space as well as in energy space. 
One thus faces a theory of two coupled fields: Qi and 
4>\. The perturbative renormalization (RG) group treat- 
ment of such a theory is well documented in the liter- 
ature^^. In the absence of electron-electron interac- 
tion (no fluctuations of </>i), the electron states in a low- 
dimensional (d = 1,2) array are localized. Accounting 
for even small fluctuations of 4>\ results in a finite phase 
coherence length. This length gets shorter at higher tern- 
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peratures, and, in the cas o 13 i 7 of d — 2, becomes of the 
order of the period of the granular array at T ~ g 2 5. 
At higher temperatures, single-particle localization is not 
important. However, as has been shown in this paper, if 
8 -C E c exp[—g], the conductivity of an array still may 
be strongly suppressed due to the formation of charge 
solitons. Those appear due to the 0-field configurations 
with non-zero winding numbers, W , or instantons. The 
latter are absent in the perturbative RG, making it in- 
adequate for the quantitative description of the insulat- 
ing phase. It is known from the single-dot models that 
the ^>-field instantons induce non-perturbative rotations 
of the Q-matrix fields. The full theory accounting for 
the single-particle localiztion effect and for the charge 
solitons, should therefore contain combined instantons 
of 4> and Q degrees of freedom™. The experience in a 
non-perturbative treatment of replica (or Keldysh) non- 
linear cr-models is relatively limited, but rapidly grow- 
ing. We should mention recent successful treatment of 
the Wigner-Dyson spectral statisticsS^i^. Even more 
encouraging is the recent realization that Anderson lo- 
calization in certain one-dimensional systems may be 
treated exactly by summing up all stationary configura- 
tions of the Q-matric field with non-zero winding num- 
ber a^?^ 4 . These developments make us optimistic that 
a quantitative treatment of the insulating phase of dis- 
ordered interacting electronic systems may be attained 
rather soon. 
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APPENDIX A: PHASE MODEL 

1. Derivation of the AES action 

We start from a prototypical model consisting of one 
dot (D) coupled to a lead (L). The system is described 
by electronic degrees of freedom and its action is given 
by 

sty] = £ SiM + StM + sM, 

where the non-interacting action Sd (Sl) of the dot 
(lead), a tunneling term St, and the charging action, S c , 



are given by, respectively, 
fP _ 

Si[i>] = / dr^(T)(5 T -/i + ^)^(r), 
Jo 

S t [tp] = J dr ^ d (t)%(t)+L(!.) , and 

ScM = E c / dr (M D (T)-q) 2 , 
Jo 

and tp(r) is an imaginary-time fermionic (Grassmann) 
field. Further, 

• Hd (Hl) is the Hamiltonian of the dot (lead). The 
notation ipiHiipi implies a sum over the internal 
Hilbert space, i.e., ipiHiipi = J d'x ^(x)ifj'0j(x). 

• T is the tunneling operator between the dot and 
the lead. Its real space representation is given by 
some matrix T(x, x') whose detailed structure we 
need not specify. 

• E c = e 2 /(2C) is the charging energy, where C 
is the self-capacitance of the dot, and Nd = 
J dx iPd(x.)iPd(x) the number operator of the dot. 

• q = V g C/e is the background charge on the dot set 
by an external gate voltage V s . 

The decoupling of the interaction part of the action is ef- 
fectuated by the introduction of a Hubbard-Stratonovich 
field V(t) with the physical significance of a voltage: 

e -S D [4>]-s c [4>] > J j)V e -S[V\-S D W,v] ^ 

where 

1 l & ft 3 
Sc[V] = — / drV 2 {T)-iq / V(t), 
c J « 

S D [lp,V] = / dT^ D {T){d T -H + H D +iV{T))i> D {T). 

Jo 

The field V can be removed from the dot action Sd by 
the gauge transformation 

^d(t) -> e l * (T) fe(r), where 4>{t) = V(t). (Al) 

However, that transformation requires some caution. 
The fermionic fields obey anti-periodic boundary con- 
ditions ipn(0) = —ip d(/3)- In order to preserve this prop- 
erty, the field </> has to fulfill the condition </>(/?) — (f)(0) — 
2irW, where W € Z. Therefore, the static contri- 
bution Vq — J dr V(t) can only be removed up to 
SV =V - 2ttTW, where W = [V /{2ttT)} is the closest 
integer to Vo/(2nT). However, in the limit of negligible 
level spacing 6 —> 0, fluctuations around 6Vq = are 
suppressed. 
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The gauge transformation (|A1|) couples to the tunnel- 
ing terms 





Since the action is quadratic, the fermionic fields can be 
easily integrated out to yield a description in terms of the 
phase field </> only. Performing the Gaussian integration 
over ip we obtain 

where the trace extends over all index spaces (time, po- 
sition, i = D,L) and G" 1 = d T — [i + Hi. 

As a next step, wc expand the ' tr In' in tunneling am- 
plitudes, 



St M = tr ( G DTe^G L f^e-^ 



(A2) 



Finally, the trace over the Hilbert spaces of dot and lead, 
respectively, leads to 



St 



1 



^K fc tr [(Ae^Ae - *)' 



(A3) 



where Kk 



i (-if 



-l^T.cS^v D v L \T a \ 2 ) k i v D {v L ) is the 
density of states of the dot (lead) and 



A(t-t') = 



sin(7rT(r - r')) 



(A4) 



(In the main text, A is given in Matsubara representa- 
tion.) Adding to this the phase representation of the 
charging action, S c [<f>] = J dr <j?/(4E c ) - iq(<f>(/3) - <l>(0)), 
we obtain the phase action of the system, S — S c + St . 

The straightforward generalization of this result to d- 
dimensional array geometries is given by Eqs. J3J), (@J. 



2. Fluctuation determinant 

In this appendix we provide details of the computation 
of the determinants resulting from the integration over 
fluctuations around instanton solutions. We will discuss 
the cases of isolated dots, one- and two-dimensional ar- 
rays in turn. 



a. Single quantum dot 

We begin by expanding the action (|T2|l to second order 
in deviations, Scj), from cj)w({z})- Since the instantons are 
saddle point configurations, there are no linear terms in 
this expansion: 

5S iast = g{5ct > \F\6<t>}, (A5) 



where the operator F is given by 

F(z) = •Aw(z) • A - i • »A w (z)A - ^Aw(z) • »A, 

and A w (z) = e i<t>w ^ Ae~^ w(z) . The dots indicate places 
for <5</>(t) in the matrix products. In a more explicit no- 
tation, 

2F(t,t') = A^'A^+A— 'KQ t 

-S(t-t') J dr, [AyJ' 1 A T1_r + A— 1 A^ T ] . 

For z = 0, the operator F is diagonal in the Matsubara 
frequency basis. Its spectrum is given by 



X (W) 



0. 



1< Irnl < \W\ 



\m\-\W\, \W\<\m\ 



For z ^ 0, the eigenbasis has a more complicated form, 
but the eigenvalues are independent of z. For W = 1, 
the basis reads 



,m+l 



1 — UZ* 



u — z 
u 



,m— 1 



1 — UZ* 



m < -2, 
m = — 1, 
m = 1, 
m > 2, 



(A6) 



where u — exp[27riTr]. [m = corresponds to a constant 
shift which is of no interest.] The quadratic action is thus 
given by 

5 inst =Go\W\- 2 m Wq + gY J ^ ) \H m \\ (A7) 



where <5</>(r) = J2 m Hm\<Pm{z))- 

Now we can evaluate the partition function taking into 
account all instanton configurations. The partition func- 
tion Z is given by Eq. JTSJl, 



,2-rriWq 



(A8) 



ir 



where Zq is the partition function in the absence of in- 
stantons. The contribution to the partition function Z\y 
corresponding to a certain winding number W consists 
of all configurations with s + W instantons and s anti- 
instanton a 28 ' 30 , where s > max{0, — W}. Here we neglect 
the weak interaction between (anti-)instantons (cf. main 
text.) Thus, 



Z\v 



E 



(2 S + wy. z 2 



s+W 



s=max{0,-W} 



(s + W)\s\ Z a 



(dz)J 



(2s+W) 



where J'^(z) is the Jacobian of the transformation to 
the collective coordinate basis {|™)} — > {z, \tp m >w)}, 



jW{z) = Adet w 



1 - Z aZ* a > 
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Furthermore, Z w is obtained by Gaussian integration 
over the massive fluctuations 8<j>, namely 



2> w 
2o 



-g oW ( Y\Z=i 9^ \ _ ( g tt m \ 



Rewriting the product Yl m (. ■ ■ ) = exp[J] m ln(. . . )], one 
finds that it is dominated by large |m|. Therefore, it is 
sufficient to keep only the first term in an expansion in 
— an upper cut-off is provided by the charging term 
in the action. Thus, 



(A9) 





( 2 E c 




- Q — e 
V T 


Finally, 




j d 2w z det M 


1 




1 - Z aZ* a > 



Putting all the components back together yields 



1 / TP P \ 2s +\W\ 



and, subsequently, 
Z 

w 
exp 



Q a 

T 



In 



El 

T 



27rg 2 ^e- g ° In ^ cos(27rg) 



(A12) 



where /,/(z) is a Bessel function. This is the result quoted 
in the main text. 



b. Fluctuation determinant of the Id array 

It is relatively straightforward to generalize the anal- 
ysis of the previous section to an array geometry. The 
main difference to the single dot case is that the fluctu- 
ation matrix F — {F k i\ has additional spatial structure, 
where the indices k, I label the dots in the array. With 
the notation 



A« = 



1 



\W k+1 
\W k+1 



w k \ = o, 

W k \ = l, 



(A13) 



where m is the Matsubara index, the fluctuation matrix 
F for the array takes the form 

F kk = A( fe - X ) + A<*>, F kk+1 = F k+lk = -A( fe ),(Al4) 

and F k i = otherwise. Here, we have put the instanton 
parameter z = 0. (This simplification is justified because 
the dominant contribution to the fluctuation determinant 
comes from high frequency fluctuations which are not 
significantly affected by the value of z.) 



It turns out to be convenient to rewrite F = 
Fq — 5F, where Fq is the fluctuation matrix on 
the flat configuration without instantons, Fq 1 = 
\m\ (24,; - 1 - 4+1,0- % contrast, SF has non- 
zero entries only at the edges of the plateau: 

SF kt = ^2 (h,i(Sk,s + 4,8+i) 

{s,s+l}edA 

— 4,«fy»+i - 4,s+i4,sJ, (Ai5) 

where dA is the boundary of the plateau, i.e., {s, s + 1} e 
dA means that \W S — W s +i\ = 1. 

Thus the fluctuation term T — detF/detfo can be 
represented as 

T = Il exp tr ln (\ m \ F ° - 6F) - tr In (\m\F { 

m 

= exp J2 tr ln f 1 - \k\ P * 1SF ) 



(A16) 



where Fq = |m|.Fo- 

The Matsubara sum starts with m = 2, i.e., there is no 
infrared divergence. The important contribution comes 
from large \m\. It is therefore sufficient to keep only the 
first term in r^r, i.e., 



N 
T ~ exp 



E 



■ tr 



F^SF 



(A17) 



Inverting the matrix Fq we find (F 1 ) k i = min{fc, 1} 



kl 
M+l 



which leads to 



tr 



F^SF 



2M 
M + l 




(A18) 



We finally note that the summation over m has to 
be cut off at large frequencies m ~ E c /gT where the 
charging energy ~ E^io^ and the dissipation action 
~ <?|a> m | become comparable. This leads to the estimate 
\ m \~ 1 ~ \ n (gE c /T). Substitution of this formula 
along with Eq. (|A18|) into Eq. (|A17|) leads to the result 

(J321- 



c. Fluctuation determinant of the 2d array 

The two-dimensional fluctuation determinant may be 
obtained by straightforward generalization of our discus- 
sion of the previous section. Concentrate for simplicity 
on the infinite system, M — ► 00, we find that (the Fourier 
representation) of the fluctuation matrix on the flat back- 
ground (no instantons) is given by Fq = \m\Fo, where 



F (p) 



2 Px 



,2 Py 



(A19) 
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To account for the presence of instantons, we again need Or, taking into account instantons and anti-instantons, 



to evaluate the quantity tr 



F^SF 



Since SF has non- 



vanishing entries only at the boundary of the island with 
W = 1, one finds 



tr 



F^SF 



= J2 {(F - l U + (F - v 

{s,s+ei}SdA 

— (F 1 ) ss +e i — (F Q )s+eii 



where {s, s + ei} e dA means that the link between sites 
s and s + ei crosses the boundary of the island, i.e., \W S — 
W s+ei | = 1. Using Eq. (|A19|I , one obtains 



tr 



F^SF 



sin 

{s.s+eijeSA p 2 



sin 



2 Py 



{s.s+eijeai 



1 _ L 

2 ~ 2"' 



where L is a number of links along the circumference 
of the island. Using the approximations l|A17(l and 
J2 m |to| _1 ~ ln(gEc/T) (which apply regardless of di- 
mensionality) we then find that the fluctuation determi- 
nant is given by Eq. (|52|l . 



Q{iu> m ) = I V V g£ + T '^ I , ) ; | e +l (A22) 



where 



s=0 v 



= Ag 2 Z w [ d 2w z jM{z) (|(^|F w |m)| 2 ) 5 . (A23) 



Furthermore, 



\(H\F w \m)t 
= ^2\i\ k '(m\ipk{z))(ipk'{z)\m) (S4>k5^-k>) S!j) 

k,k' 

and, according to the action Eq. I|A7|I . (SipkScp-k' 
Z w 5 k k'/(2g\k)- Thus, 



= 2 5 Z^ / d 2 ™z jW(z)(m|F w (z)|m). 



3. Conductance 



For small to, one finds 



To compute the conductance, we couple the action to 
source terms Ax {X = L,R), 

^ fe X=L,R L 

Physically, these terms represent the vector potential of 
an electric field coupled to the system. The Kubo conduc- 
tance, G, may be computed by taking two-fold deriva- 
tives of the partition function with respect to Ax at 
A x = 0: 



it . T 

G = lim — 3 \Q(iuj r , 

2 w^o to 



where 



1 



Q(iuj m ) = — 



S 2 Z 



Z 5A R (w m )5A L (-uj m ) 



(A20) 



(A21) 



A R =A L =Q 



Representing the partition function as a sum over wind- 
ing number sectors, one finds 



,2-rriWq 



W 



SS 



5A R (cj m ) 8A L (-cj m ) 



[n,\r ir (-:)\m) in ( l + ^ln| Zj | 2 | +Olnr). (A24) 
Finally, using that 



W\ ±A - 1—2 



(A25) 



j=i 



we obtain 



/W W -. 

d ^ z{ l + J2ln\ Zj \ 2 )l[-^r 

.7=1 .7=1 ^ 



TO / E c 

2g—Z w hrln — 
w\ \ 1 



1 + Aw 



Trlnff 



,(A26) 



where A = —ir 3 /6. 

Substituting these results into Eq. (|A22(I we obtain 
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7 °° 1 



ir 



( s + \w\y.si 



, E r 



E, 



f^e- So ln^ 

r t 



E 

2m.g(l + 2Aff 2 ^e" e ° cos(27rq)) 



2s+\W\ 



l + A{2s+\W\) 



lnff 



Recalling that w m — lixmT and substituting in 
Eq. l|A20p . one finds Eq. (|2*0j) for the conductance. 



4. Coulomb gas mapping 

We briefly review the mapping (Id Coulomb gas)^ (Id 
discrete sine-Gordon model)^. The action of the latter 
is given by 



S[6] 



El 

T 



dx 



{(<96>) 2 - 2ge- g °/ 2 cos(2tt0)},(A27) 



where we have used specific expressions for the charge 
and the fugacity matching our model parameters 
(cf. Eqs. and QQ.) 

Starting from the Coulomb gas representation, 
Eq. Ij31|l . we introduce the density variable p{x) = 
J2i <Ji8(x — Xi). Using the equality 1 = J Dp 5(p(x) — 

^2 i a i S(x-x i )) = j DO Dp exp[i9(x)(p(x) -J2i a iK x ~ 
Xi))] , we may rewrite the interaction term of the Coulomb 
gas in terms of the new variables and subsequently inte- 
grate out p. This obtains a representation in terms of the 
Lagrange multiplier field 9, 

M 



2k 

J] [ dXi e # £?,/=i I^-M- 3 <V 
•=i„ 



dx e iB ^ 



dx 



Instead of imposing strict charge neutrality, one may as- 
sume that positive and negative charges fluctuate inde- 
pendently 



£(< 

k=l 



-Q(T)/2 



2k I 



<7=±1 k=l 



(fc!) 2 



dx e l0l ~ x) 



dx e iae{x) 



dx e"*» 
(A28) 



Performing the ^-summations and relabeling 9 — > 2?t9, 
we obtain Z = J DO exp(— S[9}), where the action is 
given by $M7\ . 



APPENDIX B: CHARGE MODEL 

1. Multi— channel contacts 

In this appendix we discuss the generalization of the 
classical model derived in the main text (Sec. Ill B|) to 
N > 2 channels. For every channel a = 1 ... N of a 
multi-channel contact, a field 9 a (r) is introduced. In or- 
der to clarify the following evaluation scheme, the reflec- 
tion coefficients r are given indices specifying the direc- 
tion, contact, and channel. For simplicity, we consider 
the case of zero gate voltage q = throughout. The 
quadratic action of a d-dimensional array of M grains 
with N channels in each of the dM d contacts is then 
given by 




>,(B1) 



while the backscattering is described by 



p 

r 

S r = — ^ E E Tila / dr cos ( 27r #a,a)- 



(B2) 



Here, it is assumed that the high-energy modes E c < 
\uj m \ < D have already been integrated out. [At energies 
larger than E c , all modes are decoupled and, thus, can 
be integrated out for each channel separately] 

Only the dM d symmetric modes 0\ — J2 a @l,a couple 
to external parameters, such as gate voltages. We thus 
need to find an effective action for 0\ by integrating out 
the dM d (N — 1) asymmetric modes. To this end let us 
change variables from 0\ a (a = 1 . . . N) to 0\ and 0\ a — 

e ha - (Ol ~ £a<>a h,J)l{<* + 1) (a, a' = 1 . . .N — 1). 

While the charging term renders the symmetric fields 9\ 
massive, all the asymmetric fields 0\, a are massless. As 
a result, in the perturbation theory in powers of rn a , 
terms containing the massless fields 6\, a in the exponents 
(cosines) vanish. Rewriting the backscattering action in 
terms of the new fields, one can see that the lowest order 
non-vanishing terms are of the order Ha=i r Ha, where 
the product runs over all channels of a given contact: 
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Taking the averages 



J * — * rv 



with the actions 
a+ 1 , 



we obtain n«=i n Q <>a(^-T a ,r 2/JV for the prod- 

uct of correlators flaC ■ ■ )e, ■ The re-exponentiated ac- 

(B3) tion of thus assumes the time non-local form 



^=Ej?Efeirf+^(^ft) 2 )- f En^/ d ^ II 7t^twf cos (|Em^))}. (B4) 

1 L m ' i ce=l q'>q V " « / Q J 



As a next step, in analogy to the single-channel case 
(see Sees. II V 51 and IVB)l . one may integrate out all the 
remaining modes except the static one, 9 m =o- The pref- 



actor of the cosine-term E c J] 
malized according to 



N 



2^E c "/ is renor- 



7o 



7 (T) = 7o exp{ - ^ 5>(T a )0(T«<)><Wo} 

a, a' 

7o exp| - f(uj m )(l + -^ cosa; TO T aQ! /)|, 



a, a' >a 



where 



T 



(2M) 



7V-B q + 7r|w r „| 7r|a; m | 



(B5) 



and £ q = 4£ c £ t sin 2 (7rg i /(2M)). 

Since typical time differences T aa > — r a — T a > are of 
the order 1/T (the time integrals are dominated by the 
upper limit of integration), the cosine-term inside the 
exponent may be disregarded. Performing the Matsubara 
summation, one thus finds j(T) = 7o(T '/ 'E c ) 1 ~ 1 l d . 

Evaluating the multiple time integrations in the pref- 
actor of the cosine, one finds that the integral over the 
center-of-mass time r = ^ Q T a /N contributes a fac- 
tor 1/T, while the integration over N — 1 independent 
time differences r Q — r yields a constant cn multiplied 
by the logarithmic factor 32 \nE c /T (following simply 
from power counting.) The same logarithmic factor ap- 
pears in the framework of the phase model as a result 
of zero-mode integration. We shall not keep this log- 
arithm explicitly because all our evaluations of 7 are 
done up to a numerical factor only. As a result of these 
approximations, we reproduce the classical model with 
7 (T) ~ {T/E c y-^ d n« =1 r a = (T/E c ) 1 ~ 1 / d e~ s °/ 2 . 



2. Classical dynamics of the charge model 

Our strategy in dealing with the charge model was to 
eliminate all high-frequency degrees of freedom until only 
the zero Matsubara frequency remains. The resulting 
theory describes the classical statistical mechanics of in- 
teracting charges. This is perfectly sufficient to describe 
the thermodynamics of the system. To calculate the con- 
ductivity, however, one needs to retain the low-frequency 
(u; -c T) dynamics. The Matsubara formalism is not 
convenient for this purpose. Thus, we adopt the follow- 
ing strategy: we first switch to the Keldysh formulation 
and then integrate out the high-frequency components, 
reducing the theory to the classical, u <C T, sector only. 

For the sake of simplicity, we deal here with the 
Id single-channel model, defined by Eq. (|10l) . In the 
Keldysh formulation, the imaginary-time field 9(t) is 
substituted by the pair of real-time fields 8 cl (t) and 9 q (t) : 
originating from the symmetric and antisymmetric com- 
binations of the two branches of the Keldysh contour, re- 
spectively. The operator in the quadratic part of the ac- 
tion: 7r|oj m | — E C V 2 in Eq. (|1(J|) is transformed into a 2 x 2 
matrix in cl — q space. According to the general structure 
of the Keldysh technique^, the cl — cl element of this ma- 
trix must vanish (reflecting the fact that a field, purely 
symmetric on the two branches of the contour, has zero 
action.) The cl — q and q — cl elements are the retarded 
and advanced analytical continuations of the Matsubara 
operator: ±ttiu} — E C V 2 — > ±7r<9 t — E C V 2 . Finally, the 
q — q element is given by the fluctuation-dissipation the- 
orem (in equilibrium): 2io;coth(cj/2T). In the classical 
limit, u> <C T, it reduces to AiT. 

Integration over the high-frequency components of the 
fields leads to a renormalization of the backscattering 
amplitude, in exactly the same way as in Matsubara the- 
ory. The remaining renormalized backscattering part of 
the action, cos(27r#), transforms to cos(27r(0 + q )) — 
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cos(27r(0 d - 09)) w -47r0«sin(27r0 c/ ) + 0((6 1 ' 3 ) 3 ). An- 
ticipating that the quantum component fluctuates only 
weakly in the classical limit, we omit higher order terms 
in this expansion. As a result, the low-frequency part of 
the Keldysh action takes the form 



S = dt 



3 



[29] Ld t ef 



-4niT (0?) 



£ c V 2 0f + y2 sin(27T0f ) 



(B6) 



One can transform the last term in this action, employing 
a Hubbard-Stratonovich transformation to the partition 
function Z = J D6 cl D0« exp[iS], as 

e -^Tfdt(e]) 2 = /" I) ^ e _/dt{ a i re ?_2ie j (i)e|( i )}_ ^ 



Substituting this expression into Eq. (|B6|) . one notices 
that the resulting action depends on q -{t) only linearly. 
As a result, the integration over this field leads to a func- 



tional (5-function, imposing the following identity at ev- 
ery instance of time: 



%d t ef - E c v 2 ef 



^sin(2<) 



&(t)=0. (B8) 



This is the Langevin equation, where £,j(t) is the noise 
with the correlator that may be read out from Eq. (|B7|I , 
namely (&(*)&'(*')) = 2nTS jr 5(t - t'). 

Finally, in the multi-channel case one must renormal- 
ize the backscattering amplitude 7, as discussed in ap- 
pendix above. In addition, the coefficient in front of 
the time derivative must read n/g. This may be shown 
phenomenologically by requiring that without backscat- 
tering Eq. I|B8I) describes the dynamics of the classical 
i?C-circuit. From the bosonization perspective, this co- 
efficient is given by tt/N, where N is number of channels, 
simply because in the weak backscattering limit g ~ N . 
As a result, one recovers Eqs. (|44f) and (|45|l . used in the 
main text. 
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